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Abstract 

We derive evolution equations for the truncated Mellin moments of the parton distributions. We find 
that the equations have the same form as those for the partons themselves. The modified splitting function 
for n-th moment P'(n,x) is x n P(x), where P(x) is the well-known splitting function from the DGLAP 
equation. The obtained equations are exact for each n-th moment and for every truncation point xo £ (0; 1). 
They can be solved with use of standard methods of solving the DGLAP equations. 

This approach allows us to avoid the problem of dealing with the unphysical region x — > 0. Furthermore, 
it refers directly to the physical values - moments (rather than to the parton distributions), what enables 
one to use a wide range of deep-inelastic scattering data in terms of smaller number of parameters. 

We give an example of an application. 

PACS 12.38.Bx Perturbative calculations, 11.55.Hx Sum rules, Truncated moments 

1 Introduction 

QCD interactions between partons violate the Bjorken scaling [T|. The quark and the gluon distribution 
functions change with Q 2 according to the well-known Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
evolution equations The DGLAP equations can be solved with use of either the Mellin transform or 

the polynomial expansion in the x-space. The differentio-integral Volterra-like evolution equations change 
after the Mellin transform into simple differential and diagonal ones in the moment space and can be solved 
analytically. Then one can again obtain the i-space solutions via the inverse Mellin transform. The only 
problem is knowledge of the input parametrisation for the whole region < x < 1 what is necessary in the 
determination of the initial moments of the distribution functions. Using the polynomial expansion method, 
one deals with the parton distributions only in a limited range of the Bjorken variable: x < z < 1. This 
is very important because of the experimental constraints. The lowest value of x in present experiments is 
about 10~ 5 and the limit x — > 0, which implies that the invariant energy W 2 of the inelastic lepton-hadron 
scattering becomes infinite (W 2 = Q 2 (l/x — 1)), will never be attained. Thus the polynomial expansion 
technique allows one to avoid the problem of the unphysical region. 

An important role in DIS analyses play different sum rules, which refer to the moments of the parton 
distributions. One can compare experimental estimations for the moments of structure functions with 
theoretical predictions. Because experimental data do not cover the whole region of x, it is very useful 
to consider in the theoretical approach truncated moments of structure functions. The DGLAP evolution 
equations for the truncated moments have been discussed in |5)-|1U|. The presented methods, based on the 
approximate formulae are valid only for very small value of the truncation point (xo < 10 -2 ) and suffer from 
large errors as xo > 0.1. Besides, the not diagonal evolution equations for truncated moments, obtained in 
0-01 can be solved with a satisfactory precision for n > 2, while in a case of the first moment the precision 
is significantly worse. 

In this paper we derive the evolution equations for the truncated moments of the parton distributions, 
which are exact for every value of xo and for each n-th moment. We find that the truncated moments satisfy 
the DGLAP evolution with a modified splitting function P(x) — > x n P(x) in the Mellin convolution. This 
approach enables one to avoid uncertainties from the unmeasurable x — * region and also to use the smaller 
number of the input parameters. 

The content of this paper is as follows. In Sect. 2 we present the DGLAP evolution equations for the 
truncated moments of the parton distributions. A detailed derivation is given in the appendix. An example 
of an application to the determination of the contribution to the Bjorken sum rule is contained in Sect. 3. 
Finally, Sect. 5 contains a summary and conclusions. 
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2 Derivation of the evolution equations for the truncated 
Mellin moments of the parton densities 

The Q 2 evolution of the quark and the gluon densities is described by the well-known DGLAP equations 
0-0. Each equation of this system has the form 

^ = ^(P®s)(M). (i) 

In the above formula q(x,t) denotes the parton distribution, 

O 2 

i = ln A^— ' (2) 

^QCD 

ot a (t) is the running coupling and ® denotes the Mellin convolution: 

l 

{A9B){x) = J ( 3 ) 

X 

The splitting function P(z,t) can be expanded in the perturbative series of the a„(t). 
The n-th Mellin moment of the function f(x) is defined as 

l 

/» = / dxx n -\f(x). (4) 



Taking into account the relation 



i 



dxx n ~ L (A®B)(x) = A n B n , (5) 



one can find that the moments of the parton distributions obeys the evolution equation 

dq n (t) a s (t) _ 

-dT = "2^ -**(*)*.(*), (6) 

where the anomalous dimension y n (i) is a moment of the splitting function P(z,t) 

i 

7»(t) = J dzz"- 1 P(z,t). (7) 
o 

Eq.JBJ can be solved analytically and the parton density q(x,t) can be found via the inverse Mellin transform 

c+ioo 

<7(M) = J dnx- n q n {t). (8) 

c — ioo 

The only problem is the knowledge of the input parametrisation for the whole region < x < 1, what is 
necessary in the determination of the initial moments q„(t — to). 

Using the truncated moments approach one can avoid the uncertainties from the region x — > 0, which 
will never be attained experimentally. If we truncate the Mellin moment at xo, we obtain 



fn{x ) = J dxx n ~ x f(x). (9) 

The evolution equation Q implies for the truncated moment of q the following formula: 

i 

dq n (x ,t) _a s (t) f dxx u-i {P(g)q){Xyt)j (10) 



dt 2tt 
where 

l 



q n (xo,t) = / dxx n q(x,t). (11) 



We derive in the appendix an interesting relation for the truncated moments, namely 

1 

J dxx n - 1 (P®q)(x) = (P'®q rl )(xo), (12) 

XQ 

with 

P'{n,z) = z n P(z) (13) 

Hence ea. l|l()fl can be written as 

dq n (x ,t) a s(t) , , . , , 

J t = (P ® qn)(XO, t). (14) 

In this way we have obtained the DGLAP evolution equation for the truncated moments q n (xo,t), very 
similar to the original equation £Q| for the partons themselves. In our approach, P'(n,z) from l|131 plays a 
role of the splitting function for truncated moments. 

Now we can adapt to eg. 1141 1 known methods of solving the DGLAP evolution equation e.g. |11 | .|12 | . 
Dealing with the truncated moments allows us to study their evolution without making any assumption on 
the small- a; behaviour of the parton densities. One needs to know only the truncated moments of the parton 
distributions at the initial scale to (e.g. from the experimental data) , what constrains a number of the input 
parameters. In our approach we can also test different parton parametrisations comparing the theoretical 
predictions for sum rules, which involve q n (xo,t), with the experimental data. 



3 An example of an application 

Eq. 1141 can be rewritten in the full form 



i 



dq n (xo,t) _ a s (t) j dz , ( ^ x 



dt 2iv J z 

where 



*(■■?■«)- (?)"'("■')■ < w > 

In order to solve ea. 11511 . we use the Chebyshev polynomial expansion of P' and q n . The Chebyshev 
polynomials technique |13| was successfully used by J.Kwiecihski in many QCD treatments e.g. |14 | -|17 | . 
Using this method one obtains the system of linear differential equations instead of the original integro- 
differential ones. The Chebyshev expansion provides a robust method of discretising a continuous problem. 
More detailed description of the Chebyshev polynomials technique is given e.g. in the appendix of |10|. The 
Chebyshev polynomial approach in the case of the truncated moments is more effective than in the case 
of the parton densities themselves. This is because n-th moments for n > 1 are usually nonsingular, when 
xo — > 0. Thus we can obtain reliable results for every value of n > 1, xo and also Q 2 , independently of the 
input parton distributions. 

Here we would like to present solutions of eg. 1141 for the nonsinglet polarised structure function g\. In 
this way we are able to determine the contribution to the Bjorken sum rule |18 | .|19 | . For simplicity we use 
LO approximation. In our case 

q{x,t) = g? s {x,t) = gf{x,t) - gi(x,t), (17) 



q„(xo,t) = / dxx n g 1 (x,t). (18) 

The results for q n (xo,t) 1181 are shown in Figs. 1-3. In Fig.l we plot the truncated moments as a function 
of xo- We use two different inputs of gi S (x,to) at the low scale Qq = 1 GeV 2 : 

3 f s (Mo) ~ (l-^) 3 , (19) 
g? s (x,t ) ~ x- os (l - a:) 3 . (20) 

The evolution scale Q 2 — 10 GeV 2 . Figs. 2, 3 show the Q 2 dependence of the 1181 at xo = 0.01 and 
xo = 0.1, respectively. Here we consider again the both parametrisations 1191 . 1201 . We can also determine 
the contribution to the Bjorken sum rule 

M B SR{x\,X2,t) = J dxgi S (x,t). (21) 

For xi = 0.003, x 2 = 0.7 and Q 2 = 10 GeV 2 we obtain 6A/ SS fl = 1-23 and 1.07 for the inputs and CT . 
respectively. These results can be compared to the SMC data HJ, which yield 6A/ SS fl(0.003, 0.7, 10) = 
1.20 ±0.24 ±0.15. 
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Figure 1: Truncated n-th moments of the spin structure function g^ s as a function of the truncation point xo- 
The results are shown for the flat (|19|l (red solid) and for the singular (|20|1 (blue dotted) parametrisation. The 
evolution scale Q 2 = 10 GeV 2 
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Figure 2: Q 2 evolution of the truncated n-th moments of g^ s at xo — 0.01. The results are shown for the flat 
<|19D (red solid) and for the singular (|20|l (blue dotted) parametrisation. 
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Figure 3: Q 2 evolution of the truncated n-th moments of g^ s at xq — 0.1. The results are shown for the flat 
1)19(1 (red solid) and for the singular (|20|) (blue dotted) parametrisation. 

4 Summary 

In this paper we have derived DGLAP evolution equations for the truncated Mellin moments of the parton 
distributions. We have found that the equations closely resemble those for the partons themselves. The role 
of the splitting function for n-th moment plays P'(n,x) = x n P(x), where P(x) is the well-known splitting 
function from the DGLAP equation for the partons. The presented approach has an advantage that it can be 
successfully used for each n-th moment and for every truncation point xo £ (0; 1). The obtained equations 
are exact and can be solved with use of standard methods of solving the DGLAP equations. The polynomial 
expansion technique is in the case of the truncated moments more effective than in the case of the parton 
densities themselves. This is because n-th moments for n > 1 are usually nonsingular functions of x, when 
x — > 0. In this way one can obtain reliable results independently of the input parton distributions. 

The truncated moments approach refers directly to the physical values - moments (rather than to the par- 
ton distributions), what enables one to use a wide range of deep-inelastic scattering data in terms of smaller 
number of parameters. In this way, no assumptions on the shape of parton distributions are needed. Using 
the evolution equations for the truncated moments one can also avoid uncertainties from the unmeasurable 
very small x — * region. 

An analysis of the Q 2 QCD evolution for the truncated moments of the parton densities can be a valuable 
tool e.g. in the determination of the contribution to different sum rules from the experimentally accessible 
region. Particularly important is knowledge on the gluon contribution to the spin of the nucleon. Thus the 
truncated moments approach can be useful and interesting both from the theoretical and the experimental 
point of view. 

A Detailed derivation of the evolution equations for the 
truncated Mellin moments of the parton densities 

Here we show how to obtain eg. (1121 . The left-hand side of this relation can be rewritten in the full form 



l i 



l.h.s. = 




XQ X 



Using the Heaviside function 




we can change the order of integration in 1221 . namely 

1 1 

l.h.s.= J ^g(z) / dxx n - 1 Q(z-x)p(~y (24) 

:c icq 



Now we introduce instead of a: a new variable u, defined by 



x 

what implies 



X ° Z tr,r\ 
U = , (25) 



'■*■■■= J M J % (?)>(?)■ <*> 

X Q XQ 

As before, we use the Q function and get 

i i 

l.h.s. = [dz q(z) f ^ Q(z - u) (^) " P (^) . (27) 



z a: 

This allows us to change again the order of integration and obtain the final result 

l i 
l.h.s. = f (E2.Y p(El) I ' dzz n ^q{z)^{P' ®q n ){x )=r.h.s., (28) 



where P' is given by 1131 . 
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